Resonant inelastic x-ray scattering in one-dimensional copper oxides 
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The Cu K-edge resonant inelastic x-ray scattering (RIXS) spectrum in one-dimensional insulating 
cuprates is theoretically examined by using the exact diagonalization technique for the extended 
one-dimensional Hubbard model with nearest-neighbor Coulomb interaction. We find the following 
characteristic features that can be detectable by RIXS experiments: (i) The spectrum with large 
momentum transfer indicates the formation of excitons, i.e. bound states of holon and doublon. (ii) 
The spectrum with small momentum transfer depends on the incident photon energy. We propose 
that the RIXS provides a unique opportunity to study the upper Hubbard band in one-dimensional 
cuprates. 

PACS numbers: 78.20.Bh, 78.70.Ck, 78.66.Nk, 71.10.Fd 
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Insulating copncr oxides, when mapped onto a single- 
band component ,EJ have their occupied lower Hubbard 
band (LHB) and the unoccupied upper Hubbard band 
(UHB) separated by a Mott gap. Angle-resolved pho- 
toemission spectroscopy (ARPES) has succeeded in ob- 
serving the momentum dependence of the LHB for two- 
dimensional (2D) (Ref. 2) and ID (Ref. 3) cuprates. In 
contrast to the LHB, the UHB is far from our under- 
standing. 

The information about the UHB can be extracted from 
excitations from the LHB to UHB across the Mott gap. 
A representative technique for this is electron energy-loss 
spectroscopy (EELS), which has been apphed to both 2D 
(Ref. 4) and ID (Ref. 5) insulating cuprates. Another 
technique is resonant inelastic x-ray scattering (RIXS). 
When the incident photon energy uji is tuned through 
the Cu K absorption edge, the momentum dependence 
of the excitations from the LHB to UHB has been ob- 
served for insulating 2D cuprates, Sr2Cu02Cl2 (Ref. 6) 
and Ca2Cu02Cl2 (Ref. 7). On the theoretical side, the 
present authors predicted the dependence based on a 2D 
Hubbard model with second- and third-neighbor hopping 
terms Jj 

Typical examples of ID cuprates are Sr2Cu03 and 
SrCu02. Their magnetic susceptibilities are explain: 
by the ID Heisenberg model.El ARPES measurementaJ' 
showed that a photodoped hole cannot exist as a quasi- 
particle, but changes into two collective excitations, a 
spinon and a holon, as expected from the ID single-band 
Hubbard and t-J models. From EELS experiments, u 
it has been pointed out that an extended version of 
the Hubbard model with a moderate nearest-neighbor 
Coulomb interaction V is necessary to understand exci- 
tations across the Mott gap. It is interesting to know 
what we obtain in the RIXS spectrum prior to the ex- 
periments. 

In this study, the momentum dependence of the RIXS 
spectrum in ID cuprates is demonstrated by using the 
exact diagonalization technique for the extended Hub- 
bard model. We find the following characteristic features 
that can be detectable by RIXS experiments: The spec- 
tral weight concentrates on the narrow energy region due 



to the formation of excitons. This is also seen in the dy- 
namical charge-correlation function, although the depen- 
dence of the intensity on the momentum transfer is dif- 
ferent. We also find that the spectral weight distribution 
for a small momentum transfer depends on the incident 
photon energy. This is associated with the intermediate 
state that is characteristic of the RIXS process. 

A minimal model that can describe the LHB and UHB 
is the extended Hubbard model as used in the analyses 
of the EELS spectra.El Hereafter, we use the term "Srf- 
electron system" to stand for the LHB and UHB. The 
Hamiltonian is given by 

H3d - -t^(4,d,+i,, + H.C.) + C/^<^<^ 

i,cr i 

+vY,nfnf^„ (1) 
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where d| ^ is the creation operator of 3d electron with 



spin a at site i, nf ^—d\ ^di^„, and nf—nf-^- 



-nl^. The on- 



site Coulomb energy U corresponds to the charge-transfer 
energy of cuprates 

In the intermediate states of the Cu K-edge RIXS pro- 
cess, 3d electrons interact with a Is-core hole created by 
the dipole transition of an Is electron to 4p orbital due 
to absorption of an incident photon with energy and 
momentum K,-. This interaction is written as 



His-zd — ~Vc 
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where nf^. is the number operator of the Is-core hole 
with spin a at site i, and Vc is taken to be positive. This 
interaction Vc causes excitations of id electrons across 
the gap. The photoexcited 4p electron is assumed to go 
into the bottom of the 4p band with momentum kg and 
not to interact with either the 3<i electrons or Is-core hole 
due to the delocalized nature of the 4p orbital.Q In the 
final state, the 4p electron goes back to the Is orbital, 
emitting a photon with energy Uf aad momentmn Kj. 
The RIXS spectrum is then given byl 
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FIG. 1. Resonant inelastic x-ray scattering spectra for Cu 
i^-edge in the ID half-filled Hubbard model. The spectra of 
the elastic scattering process at AK—Q are not shown. The 
parameters used are U/t=10, V/t=0, Vc/t— 15, and T/t—l. 
The S functions (the vertical thin solid lines) are convoluted 
with a Lorentzian broadening of 0.2t. The inset is the Cu Is 
absorption spectra with a broadening of l.Ot, and the incident 
photon energy uji is set to the value denoted by the arrow. 
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where H^Hsd+His-sd, AK=K,, -K/, Aiu=uji-ujf, ^^^^ 
(Pk cr) creation operator of the Is-core hole (4p elec- 

tron) with momentum k and spin cr, |0) is the ground 
state of the half- filled system with energy Eq, \a) is the 
final state of the RIXS process with energy Ea, F is the 
inverse of the relaxation time in the intermediate state, 
and eis-4p is the energy difference between the Is level 
and the bottom of the Apz band. The momentum com- 
ponent parallel to the ID Cu chain is represented by AK 
hereafter. The values of the parameters are set to he. 
U/t=lO, Vc/t=15, and r/t=l as for the 2D cuprates.H 
The RIXS spectrum in Eq. (^) is calculated on a 14-site 
ring by using a modified version of the conjugate-gradient 
method together with the Lanczos technique. 

Figure^ shows the momentum dependence of the RIXS 
spectrum for the case of V/t—Q. The inset shows the attj 
sorption spectrum of the Is electron to the 4p orbital.B 
The incident photon energy Wi that appears in Eq. (|^) is 
set to be the value denoted by the arrow. This means that 
the states near the arrow satisfy the resonance conditions 
of Eq. (^). For small AK, the spectra spread over about 
8t. The energy region, however, shrinks with increasing 
AK, and the spectral weight concentrates on a narrow 




FIG. 2. Dynamical charge correlation function N{q, u) 
in the ID half-filled Hubbard model. The parameters used 
are U/t=10, V/t=Q (broken lines), and V/t—1.5 (solid lines). 
The (5-functions are convoluted with Lorentzian broadening 
of 0.2t. 



energy region at AK—ir. This momentum dependence is 
similar to that of the dynamical charge-correlation func- 
tion N{q,uj), shown by dashed lines in Fig. ^, which is 
explained by the nariicle-hole model with only a charge 
degree of freedom.|jO In the particle-hole model, holon 
and doublon bands separated by U describe the LHB 
and UHB, respectively. We note that the spin degree of 
freedom is less effectisia on N{q,uj) because of the wave- 
function factorizationll3 in the large-C/ limit. At the zone 
center of the interband transitions, one obtains a contin- 
uum with a width of St above the gap of U—At. The width 
of the continuum decreases with an increase of the mo- 
mentum transfer, and finally becomes an excitation with 
energy U at the zone boundary. This momentum depen- 
dence is qualitatively consistent with that shown in Fig. |l| 
and also with N{q, cu) shown in Fig. ||. It is interesting 
to note that the integrated weight of I{AK,Au!), with 
respect to Alo, is almost independent of AK in contrast 
to N(q, u) where the integrated weight is proportional to 
sin^(g/2) for the large-C/ region.li2l 

The effect of the nearest-neighbor Coulomb interaction 
V on I{AK, Alo) is shown in Fig. |^, where V/t^l.5 is 
taken according to the analyses of the EELS data.d From 
the comparison between Figs. |l| and |^, we find that the 
shape of the spectrum changes remarkably for AK>tt/2, 
accompanied by the formation of sharp peaks together 
with the shift of the spectral weight to the lower-energy 
region. This is due to excitons that appear at thejHO- 
menta satisfying a condition that V>2tcos{AK/2)Jt3 A 
similar change is also seen in N{q,uj) (see Fig.^). For 
AK<tt/2, the intensity of the spectrum near the lower 
edge is enhanced by the presence of excitonic effiects, 
which is again similar to the feature of N{q,u;)t3 We 
note that the presence of V also affects the RIXS spec- 
trum in the 2D case, although the edge stxucture, which 
is a fingerprint of the shape of the UHB, El does not de- 
pend on V. 

The spectral weight depends on the incident photon en- 
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ergy LOi. By changing uji—eis-i^p from about —21t to — 14i, 
the distribution of intensity in the lower-momentum re- 
gion shifts to the higher-energy region, as shown in Fig. ^. 
This can be attributed to the nature of the intermedi- 
ate state around the resonance condition, denoted by 
the arrows in the insets of Figs. and |^. The peaks 
at w— eis-4p~— 21i (~ U—2Vc) in the absorption spec- 
trum are characterized by states having configurations 
with the core-hole site doubly occupied by 3d electrons 
together with an empty site of id electrons. On the other 
hand, in states near uj—eis-Ap—~l'^t (~ —Vc), the core- 
hole site is singly occupied by 3d electrons and there is 
no empty site. Since the different intermediate states are 
used depending on the incident photon energy, the RIXS 
spectrum depends on the energy. The spectrum in the 
low-energy region of Acj/t^l in Fig. ^ corresponds to 
singlet excitations of the spinon pair, because the inter- 
mediate state has singly occupied configurations. 

In summary, we have examined the momentum depen- 
dence of the RIXS spectrum in ID cuprates, which pro- 
vides a unique opportunity to study the upper Hubbard 
band, using the exact diagonalization technique for the 
extended Hubbard model. We have found the follow- 
ing characteristic features of the RIXS spectrum: The 
spectrum with large momentum transfer, A/C>7r/2, in- 
dicates the formation of excitons, i.e., bound states of 
a holon and a doublon. We have also found that the 
spectral weight distribution for small momentum trans- 
fer, AK<'k/2, depends on the incident photon energy. 
The dependence is associated with the intermediate state 
that is characteristic of the RIXS process. These findings 
will be good guides for RIXS experiments on ID insulat- 
ing cuprates that will be performed in the future. 
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